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The partially attractive character of the dipole-dipole interaction leads to phonon instability in 
dipolar condensates, which is followed by collapse in three-dimensional geometries. We show that 
the nature of this instability is fundamentally different in two-dimensional condensates, due to the 
dipole-induced stabilization of two-dimensional bright solitons. As a consequence, a transient gas 
of attractive solitons is formed, and collapse may be avoided. In the presence of an harmonic 
confinement, the instability leads to transient pattern formation followed by the creation of stable 
two-dimensional solitons. This dynamics should be observable in on-going experiments, allowing for 
the creation of stable two-dimensional solitons for the first time ever in quantum gases. 






5-H 



a 



I 

O 

o 



> 

en 

00 

rn 

o 

00 

o 



X 



Recent experiments on atoms with large magnetic mo- 
ments [l|, y, y, |j| , polar molecules [5|] , and spinor Bose- 
Einstein condensates (BECs) [q are attracting a major 
interest to the physics of dipolar gases. In these gases, 
the dipole-dipole interaction (DDT) plays a significant or 
even dominant role when compared to the short-range 
isotropic interactions, which up to now have played the 
major role in the physics of ultra cold gases. The DDI are 
long-range and anisotropic, being partially attractive. As 
a consequence, dipolar gases present a rich new physics 
for Bose gases [7|, |8|, |9[ , Fermi gases [10[ , spinor gases [ll| , 



12|, and quantum information 



strongly-correlated gases 

0. 

BEC stability becomes an issue of obvious concern in 
the presence of attractive interactions. 2D and 3D ho- 
mogeneous short-range interacting condensates with at- 
tractive interactions (s-wave scattering length a < 0) are 
unstable against collapse. The presence of a trap may 
stabilize an attractive BEC for a sufficiently small num- 
ber of particles [IJ] , whereas a sufficiently large interac- 



tion energy leads again to self-similar collapse 
lowed by intermittent implosion and pattern formation 
[16| . strikingly similar to a supernova implosion 17 1. On 
the contrary, for ID BECs the attractive interactions are 
always compensated by quantum dispersion resulting in 
the formation of stable bright solitons [l8| . 

The partially attractive character of the DDI results 
in nontrivial stability conditions for dipolar BECs when 
the DDI becomes sufficiently large compared to the short- 
range interactions. In particular, 3D homogeneous dipo- 
lar BECs are unstable due to instability at very low mo- 
menta (phonon instability) |7|. For trapped 3D dipolar 
BECs this instability leads to condensate collapse resem- 
bling the above mentioned case of short-range interact- 
ing gases, although the post-collapse dynamics is quali- 
tatively different, being characterized by quadrupole-like 
patterns as recently observed experimentally [3J. 

Due to the anisotropic character of the DDI, the gas in- 
stability is largely dependent on the trap geometry. Disc- 
shape traps surpassing a critical aspect ratio have been 



recently shown to prevent phonon-like instabilities if the 
dipoles are oriented orthogonal to the trap plane [4] . In 
that case, however, a sufficiently large dipole moment 
may lead to the rotonization of the Bogoliubov spectrum 
of the dipolar BEC [SJ], which eventually leads to local 
collapses [19], or stabilized modulated density profiles in 
sufficiently small traps [20,] . 

Phonon instability in 2D traps is however present if 
the dipoles are oriented parallel to the trap plane. The 
special long-range anisotropic character of the DDI leads 
naturally to the physically relevant question of how di- 
mensionality affects the instability dynamics, and in par- 
ticular whether 2D phonon instability is necessarily fol- 
lowed by condensate collapse. In this Letter, we show 
that indeed 2D phonon instability in dipolar BECs does 
not necessarily lead to collapse, hence differing qualita- 
tively from phonon instability in 2D and 3D short-range 
interacting BECs, and 3D dipolar ones. The absence of 
collapse is explained by the formation of stable inelastic 
2D bright solitons, whose stability results from the long- 
15| fol- range character of the DDI 2l|, [H, |23| , being unstable 



in short-range interacting BECs. As a consequence, 2D 
phonon instability is accompained by the transient for- 
mation of a gas of 2D solitons, which resemble the for- 
mation of soliton trains in unstable BECs with a < 
[18J . However, contrary to the latter case, 2D dipolar 
solitons attract each other, and scatter inelastically [2J], 
eventually fusing into larger solitons, which collapse only 
if the number of particles per soliton surpasses a critical 
value. In the last part of this Letter, we show that in 
the presence of a significant in-plane harmonic confine- 
ment, 2D phonon-like instability leads to transient pat- 
tern formation followed by the creation of stable solitons. 
Interestingly, as we discuss at the end of this Letter, the 
2D instability may be easily observed in on-going exper- 
iments, allowing for the creation of 2D stable solitons for 
the first time in quantum gases. 

In the following, we consider a BEC of particles with 
mass m and electric dipole d (the results are equally valid 
for magnetic dipoles) oriented in the z-direction by a suf- 



ficiently large external field, and that hence interact via 
a dipole-dipole potential: Vd{r) = a(f{l — 3 cos^ {6)) / r^ , 
where 9 is the angle formed by the vector joining the in- 
teracting particles and the dipole interaction. The coeffi- 
cient a can be tuned within the range —1/2 < a < 1 
by rotating the external field that orients the dipoles 
much faster than any other relevant time scale in the 
system [25| . At sufficiently low temperatures the physics 
of the dipolar BEC is provided by a non-local non-linear 
Schrodinger equation (NLSE) of the form: 
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where /dr|-0(r, i)p ^ N, g — Anfi^a/Tn, a > is the 
s-wave scattering length, m is the particle mass and V 
is the external confining potential. In the following we 
use the convenient dimensionless parameter f3 — a(fi /g, 
that characterizes the strength of DDI compared to the 
short range interaction. Note that the ratio /3 may be 
easily controlled experimentally by means of Feshbach 
resonances, as recently shown in Ref. |3J. 

The 3D homogeneous {V{r) = 0) solution of Eq.([T]) 
is given by ^(f, t) = ^/n^ex.p [-i^^dt/h], where nsd de- 
notes the 3D density, and ^3^ — {g + Vd{0))n3d is the 
chemical potential, with Vd{q) = {8TT/3)a(P[3ql/\q]'^ — 
l]/2 the Fourier transform of the DDI for dipoles ori- 
ented along the z-direction. A Bogoliubov analysis of the 
3D homogeneous BEC provides the dispersion relation 
for quasiparticles: e(g) = [Ekin{q)[Ekin{q) + Eint{q)]Y^^, 
where Ekin{q) = {f^q^p + h'^ql)/2m is the kinetic energy, 
and Eint{q) = 2(3 -I- Vd{q))n3d is the interaction energy, 
which includes both short-range and dipolar parts. Note 
that Vd{q} may become negative for some particular di- 
rections, and hence for low momenta phonon instability 
is just prevented if — 3/87r < /3 < 3/47r. If a > 0, phonons 
with q lying on the xy plane are unstable if /3 > B/Att, 
while for a < phonons with q along z are unstable if 
f3 < —3/8tt. In both cases, the dipolar BEC is unsta- 
ble against local collapses. In the presence of trapping, 
this instability leads to a global collapse of the BEC, and 
may be geometrically prevented in sufficiently pancake- 
shapped traps, as shown recently in Ref. '^j. 

In the following we show that the phonon instability 
may become crucially different in 2D dipolar condensates. 
We first consider the case of homogeneous 2D dipolar 
BECs in the xy plane, which are strongly confined by an 
harmonic potential V = (l/2m)a;^z^ along the transver- 
sal z-direction, such that the system can be considered 
"frozen" into the ground state (t)o{z) of the harmonic os- 
cillator in the z direction, and hence the BEC wave func- 
tion factorizes as ^(r) = ^p{p)(j)o{z). The dipole orien- 
tation with respect to the trap plane plays of course a 
crucial role. In the following we discuss the two extremal 



configurations, in which the dipoles are normal to the 
xy plane (_L-configuration) and parallel to it (along x) 
(||-configuration). Employing the above mentioned fac- 
torization, the convolution theorem, the Fourier trans- 
form of the dipole-potential and integrating over the z 
direction, we arrive at the 2D NLSE 2l|: 
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where k is the xy-momentum, n the Fourier transform of 
|f/'(p)P, ^2_d(k) = 2 — 3^/7rKe'^ erfc(K) (_L-configuration) , 
and /i2d(k) = -1 + 3y/^{Kl/k)e''\Tic{K) (||- 
configuration), with erfc(K) the complementary error 
function and l^ — i/^i/mw^. 

The homogeneous solution of ([2]) is i/'(P; t) — 
\Jn2d exp [— «/J.2di/^] 1 where n2d denotes the homoge- 
neous 2D density, and [iid — gn-id^X + 87r/3/3)/v27rZz 
is the chemical potential for _L-configuration and \xid = 
gn2di^ — At^ P / 3) / \/2ttI z for | |-configuration. The corre- 
sponding Bogoliubov equations provide the dispersion re- 
lation for elementary excitations of 2D BEC 22| : 
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where Ekin = h^\k\'^ /2m. 

We discuss first the _L-configuration. From ([3]) we see 
that the excitations with low momenta (fc ^ 0), become 
purely imaginary if /3 < — S/Stt, i.e. phonon instability 
demands the tuning of the DDI, being absent without 
tuning. As discussed in Ref. [2l| stable 2D dipolar bright 
solitons are possible for g > if /? < — 3/87r. In the fol- 
lowing we call these solitons isotropic solitons (since they 
are isotropic on the xy-plane) to distinguish them from 
the solitons introduced in Ref. [23| , which we discuss be- 
low. As observed, the phonon-instability and the soliton- 
stability conditions are exactly identical. This fact, far 
from being accidental, shows that the phonon-instability 
and 2D soliton formation are intrinsically linked. 

Using Eq. Q , we have numerically studied the phonon 
instability by employing periodic boundary conditions. 
In all our simulations, we check that /l2d = M2d/Sa;^ < 1 
to ensure the 2D character of the problem. Starting from 
an homogeneous solution, our real time simulations show 
the dynamical formation of a soliton gas, as shown in 
FigH] (top). The long-range character of the DDI leads 
to an attractive soliton-soliton interaction (since /3 < 0). 
As discussed in Ref. [2J| the soliton-soliton scattering is 
inelastic, and for low relative velocities (as it is the case in 
the soliton gas) soliton fusion is expected. This is indeed 
observed, and hence the number of solitons decreases 
in time, whereas the number of particles per soliton in- 
creases. As long as the problem remains 2D, the isotropic 







Figure 1: 2D soliton gas after phonon instability of an homo- 
geneous dipolar BEC. For both the figures /i — —0.2 and 
g/{V2Tvl^) = 10. (top) ^-configuration, P = -0.2, t = 
484/cJz; (bottom) ||-configuration, /? — 0.3, and t — 108/ luz- 
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Figure 2: Stability diagram of an anisotropic soliton as a 
function of /3 and gcr = gNcr/V^lz, where for A'^ > Ncr 
the soliton is unstable against collapse even for f3 > 3/47r. 



solitons remain stable (contrary to the anisotropic soli- 
tons discussed below). However, the density within a 
given sohton may eventually become large enough to vi- 
olate the 2D condition jj,2d < 1. In that case the soliton 



becomes unstable against 3D collapse [21 1. 

The I |-configuration shows a relatively similar physics. 
From Eq. ^ we see that e{k -^ 0) becomes unstable 
a (3 > 3/Att. This condition is exactly that found re- 
cently in Ref. [23\ for the stability of 2D bright solitons. 
Following Ref. 23] we call in the following these solitons 
anisotropic solitons, since the solitons are more elongated 
along the dipolar axis (in our case the x axis) . Before con- 
sidering the instability dynamics in the | |-configuration, 
we briefly comment on the properties of the anisotropic 
solitons, which differ to some extent from those of the 
isotropic solitons. As mentioned in Ref. [23| a simple 
Gaussian Ansatz on the xy-plane with unequal widths 
Lx.y, shows the appearance of a minimum in the en- 
ergy functional E{Lx,Ly) for some equlibrium soliton 



widths Lxo and LyQ. However, one may show that for 
a given (3 > 3/47r, there is a critical universal value 
QcriP) = gN^r/V^lz (scc Fig. [l]) such that for N > Ncr 
the minimum of E{Lx,Ly) dissappears. As a conse- 
quence, contrary to the case of isotropic solitons, which 
are always stable as long as the soliton is 2D, there is a 
critical number of particles per soliton, Ncr, which de- 
creases for larger /3. Beyond this number the 2D soliton 
collapses. This result is also verified by a direct simula- 
tion of the 3D Gross-Pitaevskii equation ([T|) . 

Using Eq. ^, and starting once more from an homo- 
geneous gas, we have observed the dynamical formation 
of a gas of anisotropic solitons, after initial formation of 
stripes along the dipole axis (Fig. [U bottom). Due to 
the anisotropic nature of DDI in the plane, the inelastic 
dynamics of anisotropic solitons is relatively more com- 
plicated. Solitons along the x axis attract each other, fuse 
together and become unstable against 2D collapse once 
they surpass the critical number of particles per soliton. 

In the final part of this Letter, we consider a 2D dipo- 
lar gas in the presence of an harmonic confinement on the 
xy-plane, V{p) = mLu'^p'^/2. As an initial condition for 
our time evolution, we have considered a dipolar BEC in 
the Thomas- Fermi regime (sufficiently large g) , obtained 
for the _L-configuration with /3 > (no phonon instabil- 
ity) by imaginary time evolution of Eq. ([2]). In order to 
trigger the phonon-like instability in our real time evolu- 
tions, at time i = we may either switch /3 < (which 
keeps the polar symmetry) or tilt the dipole switching 
into the | [-configuration (hence violating the polar sym- 
metry). We call these two cases, cases (i) and (ii). 

Figs. [3] illustrate the case (i). The phonon-like insta- 
bility leads to transient ring-shapped structures. The 
harmonic confinement leads to successive shrinkings and 
expansions of the rings, which oscillate in the trap. This 
process continues until the inner rings merge into a sin- 
gle excited dipolar bright soliton. We have observed that 
for sufficiently shallow traps the ring structures develop 
azimuthal instability, which leads to 2D bright solitons, 
recovering, as expected, the homogeneous case. 

The case (ii) is illustrated in Figs. [H An initial 
anisotropic ring-like structure breaks into a pair of 
anisotropic clouds, which eventually merge into the trap 
center forming a single stable bright soliton, which as ex- 
pected is anisotropic, being more elongated in the dipole 
x-direction. Note that both for the cases (i) and (ii) the 
non-dissipative character of the problem implies that the 
solitons are produced in an excited internal state. 

Summaryzing, we have shown that the phonon insta- 
bility of a 2D dipolar BEC differs qualitatively from 2D 
and 3D short-range interacting gases, and 3D dipolar 
BECs. Contrary to these cases, the phonon instability 
does not necessarily lead to the collapse of the gas, but 
on the contrary leads to a transient regime characterized 
by the formation of a gas of attractive inelastic 2D bright 
solitons, which eventually undergo fusion, leading to the 
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Figure 3: Trapped BEC in the _L-configuration. In the simula- 
tions, Ljp = O.OOIlo^, g/{V2^lz) = 200 and /3 = 0.3 {t < 0) and 
—0.3 {t > 0). (left) Formation of ring structures (i = 810/lJz); 
(right) Final soliton at the trap centre (i = 7500/llIz). 




Figure 4: Trapped BEC after tilting at f = O from _L to || 
configuration, with ujp — O.OOIljz, g/{V2Tvlz) = 200 and (3 = 
0.3. (left) Anisotropic ring formation (i — 1092/a;z); (right) 
Final formation of an anisotropic soliton (f = 1281/a;z). 



creation of a single excited stable bright soliton. If the 
dipoles are perpendicular to the trap plane these solitons 
are stable as long as the gas remains 2D, whereas if the 
dipoles are parallel to the trap plane the (anisotropic) 
solitons may become unstable even in 2D for a critical 
number of particles per soliton. Finally, we have shown 
that the instability in the presence of an harmonic con- 
finement is followed by the creation of transient ring-like 
and anisotropic patterns, which eventually lead to the 
creation of a single excited 2D soliton. 

Finally, we comment on possible experimental realiza- 
tions. Phonon instability for dipoles lying on the trap 
plane (and hence anisotropic solitons [23|) does not de- 
mand the use of tuning, which may be experimentally 
demanding [26|. Hence, 2D phonon instability should be 
relatively easy to observe in on-going experiments. A 
possible path would be to prepare a 2D stable BEC ori- 
enting the dipoles perpendicular to the trap plane, and 
then suddenly re-orient the dipoles on the plane. An 
equivalent alternative, more feasible experimentally i 26j . 
would be to maintain at any time the dipole on the trap 
plane, but drive the instability by increasing /3 using Fes- 
hbach resonances as in Ref. 3]. As we have shown, the 
2D instability would allow for the creation of 2D bright 
solitons for the first time ever in quantum gases, since 
they are fundamentally impossible in absence of DDL 
2D solitons would offer new scenarios for the matter-wave 



soliton scattering, as e.g. 2D inelastic scattering and spi- 
raling, similar to that in photorefractive materials [271] . 
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